Mott transition in multi-orbital systems 



A. Liebsch 

Institut fur Festkorperforschung, Forschungszentrum Julich, 52425 Julich, Germany 

(Dated: February 2, 2008) 



m 
o 
o 

(N 

c 

(N 



i 

CO 

-4-* 



I 

C 

o 
o 



> 

(N 

O 
m 
O 

c3 



c 

o 
o 



X 
S3 



Metal insulator transitions driven by local 
Coulomb interactions are among the most fasci- 
nating phenomena in condensed matter physics 1 . 
They occur in a large variety of transition metal 
compounds 2 . Most of these strongly correlated 
materials consist of valence bands derived from 
electronic d shells where intra- and inter-orbital 
Coulomb interactions are equally important and 
where the crystal structure splits the valence 
bands into narrow and wide subbands. A funda- 
mental question is whether these systems exhibit 
a common Mott transition, implying all subbands 
to be either metallic or insulating, or successive 
orbital dependent transitions, implying a coexis- 
tence region with metallic and insulating behavior 
present in different subbands. Using the dynami- 
cal mean field theory 3 we show that inter-orbital 
Coulomb interactions lead to a single Mott tran- 
sition. Nevertheless, the subbands exhibit more 
or less strongly correlated excitation spectra in 
the metallic phase and different band gaps in the 
insulating phase. 

A hallmark of strongly correlated materials is their ge- 
ometrical, electronic and magnetic complexity. The va- 
lence bands typically involve oriented electronic orbitals, 
giving rise to a variety of highly anisotropic properties. 
For instance, the layer perovskite S^RuCU consists of 
a wide, two-dimensional d xy band coexisting with nar- 
row, nearly one-dimensional d xz , yz bands*. This sys- 
tem exhibits unconventional p-wave superconductivity 5 , 
but iso-electronic replacement of Sr by Ca induces 
a Mott transition to an antiferromagnetic insulator 6 . 
Similar splittings into different subbands occur in the 
classic Mott insulators V0 2 ^ and V2O3MI in both 
the low-temperature insulating (monoclinic) and high- 
temperature metallic (rutile, corundum) phases, in lay- 
ered organic superconductors^, fuUcrcncs 1 --, and many 
other compounds^. 

The nature of the paramagnetic metal insulator tran- 
sition in a multi-band environment involving narrow and 
wide subbands is not yet understood; in particular, it 
is not clear whether all bands undergo a common tran- 
sition at the same critical Coulomb energy^ or whether 
different subbands generate transitions at successive crit- 
ical energies^*. In the first scenario narrow subbands 
force wide subbands to become more strongly correlated 
and the more pronounced metallicity of wider subbands 
makes narrow bands less correlated than in the isolated 
case. The critical Coulomb energy U c of the interacting 
system lies between the critical U C i of the isolated sub- 
bands. Coexistence of metallic and insulating behavior in 



different subbands does not occur. In the second scenario 
the different relative importance of Coulomb correlations 
in subbands of different widths gives rise to separate Mott 
transitions and to coexistent metallic and insulating be- 
havior in different subbands. In this case the one-electron 
properties imposed by the different band widths remain a 
determining factor for the system's excitation spectrum 
and are not entirely superseded by correlations. 

While local Coulomb interactions in multi-band sys- 
tems have been investigated previouslyi£*i&, most stud- 
ies treated the special 'isotropic' case of identical orbitals 
where the question of one or several Mott transitions does 
not arise. Here we consider the opposite 'non-isotropic' 
situation where the multi-band system consists of differ- 
ent subbandsii. In the absence of inter-orbital Coulomb 
interactions these bands would exhibit Mott transitions 
at different critical energies f7 c j. As shown below, the 
striking consequence of inter-orbital interactions is to en- 
force one common metal insulator transition for all sub- 
bands at the same critical U c . Remarkably, however, in 
the joint metallic region different subbands reveal exci- 
tation spectra with varying correlation signatures. Sim- 
ilarly, the insulating phase exhibits a 'multi-gap' spec- 
trum, in close analogy to the a- and 7r-type supercon- 
ducting gaps observed in the quasi-two-dimensional com- 
pound MgB^. 

Let us consider the paramagnetic metal insulator tran- 
sition in a two-band Hubbard model with narrow and 
wide subbands. These might represent, for instance, the 
t2 g orbitals of a layer-perovskite material or of VO2. For 
illustrative purposes these bands are assumed to be half- 
filled and to have elliptical densities of states Pi(ui) of 
widths W\ = 2eV and W2 = 4eV. Single-band systems 
of this type have been investigated extensively in the 
past and their metal insulator transitions are rather well 
understood^. To account for Coulomb correlations we 
use the dynamical mean field theory (DMFT) in com- 
bination with the multi-orbital Quantum Monte Carlo 
(QMC) method 3 -. The quantum impurity problem is 
solved at finite temperature T for various intra- and inter- 
orbital Coulomb energies U and U' = U — 2 J, where J 
is the Hund's rule exchange integral. 

Two criteria are employed to determine the stability of 
the Fermi liquid state as a function of Coulomb energy. 
First, we calculate the subband quasi-particle weights 
Zi = 1/(1 — 9Rc £j(cj)/<9cj| w= o)j where the derivative 
of the real part of the self-energy component £i(u>) is ap- 
proximated by ImEj(zcjo)/^o with — tt//3 — itIcbT 
the first Matsubara frequency. The second criterion 
is obtained from the imaginary-time Green's functions 
at r = P/2: Gi((3/2) = JdujF(uj) ImGi(u>)/n, where 
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FIG. 1: (a) Quasi-particle weights Zi and (b) normalized 
imaginary-time Green's functions Gi(/3/2) of two-band sys- 
tems as a function of U for T — 125 meV. Solid symbols: 
isotropic Wi,Wi and W%, W2 systems; open symbols: non- 
isotropic W\,W2 system. Dotted lines: extrapolation of Zi 
and Gi(f3/2) to estimate critical Coulomb energies. 

F(uj) = 0.5/cosh(/3u;/2) is a distribution of width w = 
41n(2 + y/S)/P centered about the Fermi level. While 
Zi specifies the quasi-particle weight of the i th subband 
at Ep, Gi(P/2) represents the integrated spectral weight 
within a few ksT of Ep, i.e., it includes low- lying exci- 
tations. Since ^i(iui n ) and Gi(r) are derived directly 
from the QMC calculation Zi and Gi(f3/2) are avail- 
able without having to evaluate the real-frequency spec- 
tral distributions ImGj(w). Below we present normal- 
ized quantities G l ((3/2) = G 4 (/3/2)/Gf =°(/3/2) so that 
Zi = Gi([3/2) = 1 in the non-interacting limit. 

Figure 1 shows the variation of Zi and Gj(/3/2) as a 
function of U for T = 125 meV. The exchange integral 
is J = 0.2 eV—. Evidently the quasi-particle weights 
near Ep diminish as a result of Coulomb correlations. 
At a given U, however, Z\ and Gi in the non-isotropic 
W\ , W2 system are larger than in the degenerate Wi , Wi 
case. Similarly, Z 2 and Gi in the W\,Wi system are 
smaller than in the W2, Wi case. Thus, as a consequence 
of inter-orbital interactions the narrow (wide) subband is 
less (more) correlated than in the isotropic case. 

Since Gi((3/2) represents the spectral weight within 
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FIG. 2: Quasi-particle weights Zi of two-band systems as a 
function of U, calculated within iterated perturbation theory 
at T — 0. Solid curves: isotropic Wi, Wi and W2, W2 systems; 
dashed curves: non-isotropic Wi,W2 system. Dotted lines: 
extrapolation of Zi to estimate critical Coulomb energies. 



several ksT of Ep its reduction at small U is weaker 
than that of Zi, while in the critical regions its decay is 
more abrupt. Close to the metal insulator transitions Zi 
and Gi show a rounding off caused by the finite temper- 
ature and by critical slowing down. Within this slight 
uncertainty Zi and Gi yield consistent critical Coulomb 
energies: U c \ w 3.0 eV for the narrow two-band sys- 
tem, U C 2 ~ 6.1 eV for the wide two-band system, and 
U c w 4.8 eV for the mixed systemSi. 

The key point of these results is the fact that the 
system involving narrow and wide subbands exhibits a 
single Mott transition at an intermediate U c such that 
U c i < U c < U C 2- At small U correlations reduce the 
quasi-particle weight in the narrow subband much more 
rapidly than in the wide band. For U > U c /2, however, 
Z\ and Zi begin to converge again and decrease to zero 
spectral weight at the same U c . The same behavior is 
shown by the integrated weights Gi((3/2). There is no 
evidence for a coexistence of insulating and metallic be- 
havior in different subbands 22 . On the other hand, since 
Z\ < Zi and G\ < Gi for U < U c , the narrow band 
is always more strongly correlated than the wider coun- 
terpart. Analogously, the insulating phase exhibits two 
excitation gaps (see below). 

Multi-band QMC calculations at lower temperatures 
(T = 62meV and 31 meV) confirm the picture discussed 
above. To verify that the same conclusion holds also 
in the low temperature limit we have performed anal- 
ogous two-band calculations at T = within the self- 
consistent iterated perturbation theorji 2 ^ (IPT), where 
the self-energy is calculated at real frequencies within 
second-order perturbation theory. The results for the 
quasi-particle weights Zi are shown in Fig. 2. They are 
consistent with the scenario for T > 0: The isotropic 
two-band systems have widely different critical Coulomb 
energies: U c \ w 2.3 eV and U c i w 4.3 eV. In contrast, 
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FIG. 3: Quasi-particle spectra of non-isotropic two-band sys- 
tem for T = 125 meV (Ef = 0): (a) metallic phase, (b) insu- 
lating phase. The solid and dashed curves denote the spectra 
of the narrow and wide subbands, respectively. The bare den- 
sities of states are represented by the dot-dashed curves. 

the non-isotropic system undergoes a Mott transition at 
U c w 3.2 eV. The specific values of these Coulomb ener- 
gies differ from those in Fig. 1 since the IPT amounts to a 
more approximate treatment of correlations, and because 



of the different temperatures. 

Although inter-orbital Coulomb interactions lead to a 
common metal insulator transition in the narrow and 
wide subbands, their excitation spectra differ qualita- 
tively both in the metallic and insulating phases. This is 
illustrated in Fig. 3 where we plot quasi-particle distribu- 
tions derived using the Maximum Entropy method—. For 
U < U c both bands are metallic. The narrow band spec- 
trum Ni(lo), however, is more strongly correlated than 
N2(ui) for the wider band: The spectral weight at Ef 
is much reduced and the Hubbard peaks are more pro- 
nounced. Similarly, in the insulating phase the excitation 
gap of the narrow band is larger than that of the wide 
band. Note that the W\ , W\ system at U — 4 eV is al- 
ready insulating, while the W2,W% system is much less 
correlated than the A^w) spectrum shown in Fig. 3(a). 
Conversely, at U = 6 eV the W\ , W\ system is far in the 
insulating range, whereas the W2, W2 system is close to 
its metal insulator transition. 

These predictions can be tested using angle resolved 
photoemission spectroscopy. In the metallic phase of lay- 
ered perovskites, for example, the d xy and d XZjVZ sub- 
bands should exhibit coherent peaks of different spectral 
weights as well as different incoherent satellite features. 
At U c the coherent peaks should vanish simultaneously. 
In the insulating phase, the subbands should exhibit dif- 
ferent band gaps and Hubbard bands. 

For illustrative purposes we have considered the special 
case of half-filled, symmetric subbands. The consistency 
of the picture discussed here with the results obtained 
for Ca2- 2; Sr :I .Ru0^i^ suggests that the same conclusion 
holds for non-symmetric multi-band systems. In partic- 
ular, it should be valid for the Mott insulators VO2 and 
V2O3 and for other transition metal compounds. 

Acknowledgements 

I am grateful to A. Bringer for discussions and to O. 
Gunnarsson and D. Vollhardt for comments. I also like 
to thank A. Lichtenstein for the DMFT-QMC code. 

Electronic address: a.liebsch@fz-juelich.de 



1 Mott, N.F. Metal Insulator Transitions (Taylor and Fran- 
cis, London, 1990). 

2 Imada, M., Fujimori, A. & Tokura, Y. Metal insulator 
transitions. Rev. Mod. Phys. 70, 1039-1263 (1998). 

3 For a review, see: Georges, A., Kotliar, G., Krauth, W. & 
Rozenberg, M.J. Dynamical mean-field theory of strongly 
correlated fermion systems and the limit of infinite dimen- 
sion. Rev. Mod. Phys. 68, 13-125 (1996). 

4 Oguchi, T. Electronic band structure of the superconduc- 
tor Sr 2 Ru0 4 . Phys. Rev. B 51, 1385-1388 (1995). 

5 Maeno, Y., Rice, T.M. & Sigrist, M. The Intriguing Super- 
conductivity of Strontium Ruthenate. Phys. Today 54(1), 
42-47 (2001). 

6 Nakatsuji, S. & Maeno, Y. Quasi-Two-Dimensional Mott 
Transition System Ca2-a:Sr x Ru04. Phys. Rev. Lett. 84, 
2666-2669 (2000). 



7 Goodenough, J.B. Direct Cation-Cation Interactions in 
Several Oxides. Phys. Rev. 117, 1442-1451 (1960). 

8 Zylbersztejn, A. & Mott, N.F. Metal-insulator transition 
in vanadium dioxide. Phys. Rev. B 11, 4384-4395 (1975). 

9 Rozenberg, M. J., Kotliar, G., Kajueter H., Thomas, G. 
A., Rapkine, D. H., Honig, J. M. & Metcalf, P. Optical 
Conductivity in Mott-Hubbard Systems. Phys. Rev. Lett. 
75, 105-108 (1995). 

10 Ezhov, S. Yu., Anisimov, V. I., Khomskii, D. I. & 
Sawatzky, G. A. Orbital Occupation, Local Spin, and Ex- 
change Interactions in V2O3. Phys. Rev. Lett. 83, 4136- 
4139 (1999). 

11 Lefebvre, S. et al. Mott Transition, Antiferromagnetism, 
and Unconventional Superconductivity in Layered Organic 
Superconductors. Phys. Rev. Lett. 85, 5420-5423 (2000). 

12 Takenobu, T., Muro, T., Iwasa, Y. & Mitani, T. Anti- 



4 



ferromagnetism and phase diagram in ammoniated alkali 
fulleride salts. Phys. Rev. Lett. 85, 381-384 (2000). 
This was anticipated by preliminary results obtained in: 
Liebsch, A. Absence of orbital dependent Mott transition 
in Caa-^Sr^RuCU. Europhys. Lett. 62, xxx (2003). The 
results were, however, not fully conclusive due to uncer- 
tainties caused by the Maximum Entropy method. 
Anisimov, V.L, Nekrasov, I.A., Kondakov, D.E., Rice, 
T.M. & Sigrist, M. Orbital-selective Mott transition in 
Ca2-*Sr x Ru0 4 . Euro. Phys. J. B 25, 191-201 (2002). 
See: Gunnarsson, O., Koch, E. & Martin, R. M. Mott 
transition in degenerate Hubbard models: Application to 
doped fullerenes. Phys. Rev. B 54, R11026-11029 (1996); 
and: Ono, Y., Potthoff, M. & Bulla, R. Mott transitions in 
correlated electron systems with orbital degrees of freedom. 
Phys. Rev. B 67, 035119-1-6 (2003) and references herein. 
Florens, S., Georges, A., Kotliar, G. & Parcollet, O. Mott 
transition at large orbital degeneracy: Dynamical mean- 
field theory. Phys. Rev. B 66, 205102-1-13 (2002). 
Liebsch, A. & Lichtenstein, A. Photoemission Quasiparti- 
cle Spectra of Sr 2 Ru0 4 . Phys. Rev. Lett. 84, 1591-1594 



(2000). 

18 Suoma, S. et al. The origin of multiple superconducting 
gaps in MgB 2 . Nature 423, 65-67 (2003). 

19 See: Bulla, R., Costi, T. A. & Vollhardt, D. Finite- 
temperature numerical renormalization group study of the 
Mott transition. Phys. Rev. B 64, 045103-1-9 (2001) and 
references herein. 

20 Below U — 1 eV, J is linearly reduced to zero. 

21 Ud and U C 2 do not differ exactly by a factor of two since 
/3 and J are kept constant. 

22 As shown in Ref. 16 for the isotropic case, at very low T 
the insulating phase begins at a U' c slightly below the dis- 
appearance of the metallic phase. Similar behavior might 
occur for non-isotropic multi-orbital systems. 

23 Georges, A. & Kotliar, G. Hubbard model in infinite di- 
mensions. Phys. Rev. B 45, 6479-6483 (1992). 

24 Jarrell, M. & Gubernatis, J. E. Bayesian inference and the 
analytical continuation of imaginary-time quantum Monte 
Carlo data. Phys. Rep. 269, 133-195 (1996). 



